Spin dynamics in the XY model 
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Abstract 

We study the evolution of entanglement, quantum correlation and classical correlation for the 
one dimensional XY model in external transverse magnetic field. The system is initialized in the 
full polarized state along the z axis, after annealing, different sites will become entangled. We 
study the three kinds of correlation for both the nearest and the next-nearest neighbor sites. We 
find that for large anisotropy parameter the quantum phase transition can be indicated by the 
dynamics of classical correlation between the nearest neighbor sites. We find that the dynamics 
of entanglement for both the nearest and next-nearest neighbor sites show significantly different 
behaviors with different values of magnetic field. We also find that the evolution of quantum 
correlation and classical correlation of the nearest neighbor sites are obviously different from those 
of the next-nearest neighbor sites. 
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As an ideal model used to describing phase transition on the early days in statistical 
mechanics, the XY model has been studied intensely for almost one century. Although it 
is so simple that it can be calculated analytically, it catches the most important elements 
of some physical system, such as materials Cs(H 1 _ a .D x ) 2 P04, PbH 1 _. I .D, r P04 [1], etc. So 
it could be used to depict the actual phase transition very well. It also became the test 
bed equipped to describe more practice physical system, and to examine the performance of 
different approximation methods j^j]. With the birth of quantum information theory, it was 
found that the XY model could be used as a quantum information channel. So it became 
the focus of theoretical and experimental studying again {^J , and many quantum information 
process has been proposed with the XY model 4t6|. 

There were many works on quantum information process in the XY and the Ising model, 
such as the transfer of an unknown quantum state along the chain {4, 5, 7|. As an important 
resource in quantum information process, it was found that entanglement could be used as 



an quantity to indicate quantum phase transition 



9] . Recently, it has been reported that 



dynamics of the nearest-neighbor entanglement may be used as an indicator of quantum 



phase transition in the transverse Ising model 



a special part of the quantum correlation 



101 ]. In fact, quantum entanglement is only 
14j . there is other quantum correlation that 



is not caught by the quantum entanglement, and it is found to be important in quantum 



information task 



15 



5Hl7j]. More correlation indicates more information, so it is not trivial to 
study the evolution of the quantum correlation. 

In this paper, we study the entanglement, quantum correlation and classical correlation 
dynamics of both the nearest and next-nearest neighbor sites in the one dimensional trans- 
verse XY model with the initial state been the full polarized state along the z axis. For 
different values of anisotropy parameter, we find that the entanglement takes on distinctive 
behaviors with different values of magnetic field. We find that for large value of anisotropy 
parameter the first maximum of classical correlation between the nearest neighbor sites peaks 
around the critical point. That may be an indicator of quantum phase transition. Further 
more, we find that behaviors of the three kinds of correlation for the nearest neighbor sites 
are obviously different from those for the next-nearest neighbor sites. 

The Hamiltonian of the one dimensional XY model in transverse magnetic field reads as 



N 



N 



H = -A^l + 7)^^+1 + (1 - l)S*SV +1 ] - 



(1) 



i=l 



8=1 



where Sf (a = x, y, z) are the spin 1/2 matrixes at site i, and N is the number of sites. We 
assume periodic boundary conditions Sf = S^ +i , and the anisotropy parameter 7 divides 
the whole parameter space into two classes, for 7 = 0, it is the XX model, and < 7 < 1 
belongs to the Ising universality class. Here we will only concentrate on the thermodynamic 
limit N = 00, at which the Hamiltonian system undergoes a quantum phase transition at 
A c = 1. 

This model could be solved analytically by mapping the spin half system into the spinless 
fermion system via the Jordan- Wigner transformation 



k=l k=l 



k=l k=l 



= - \ ( 2 ) 

where bj (bj) is the creation (annihilation) operator of the spinless fermion. After appli- 
cation of the Fourier and Bogoliubov transformation rj = -^Y2i eM ( a A + i Piety), with 
parameters ah = Afc-(i+Acosfc) « _ 7Asinfc ^ e even ^ ua v]y Hamiltonian can 

be diagonalized as 

H = YM4vk- 1 -), (3) 

k 



with the eigenenergy spectrum A k = ^/(l + Acos/c) 2 + A 2 7 2 sin 2 /c. 

The time evolution of the spin system in Eq. (1), can be given by the evolution of the 
spinless fermion operator bj(t) [isl 

= >v7 J2^ cosk ^ 1 - j) ^ ~ ^smk(l - j)d k (t)b\], (4) 



k.i 

_iAtt 



with time dependent coefficients Ck(t) = e kt — 2z/3|sinAfct, dk(t) = ak/3kS\nAkt. 

We choose the initial state that all spins are polarized along the z axis, p(0) = |iV)(JV|, 
with \N) = I t)i---l t)jv, which is ground state of Eq. (1) when a strong magnetic field is 
applied along the z axis. Hence, there is no quantum correlation and classical correlation 
among any partition of the system. Due to the translational symmetry of the system, we 
choose site 1 as our first site when we study the correlation between different sites. In order to 
obtain the evolution of entanglement, quantum correlation and classical correlation, we need 
to get the reduced density matrix first. So we will have to calculate two-body spin correlation 



functions, it is equal to calculate the fermion correlation function (b±...)t = (&i(i)---)o- From 
Wick theorem, we can learn that the correlation function will be nonzero for only even 
number of creation and annihilation operators. 

For the first two spins the reduced density matrix in the basis of | "ft), | tl), I It)) I l-l) > 
is given by 



[p\i o o p\{\ 

o p\\ pE o 

o pII pU o 

\p\i o o p\lj 



(5) 



where [1 



p\l = (& t 1 (t)6l(t)4W&2(t))0=(& t l(t)&l(t))0(4W&2(t))0 

- (#)#)}o<&i(t)& 2 (*))o + (b\(t)b 2 (t)) 
x (h(t)bi(t)) , 

pII = pII = (b\(tMt)b 2 (t)bl(t)) = <&t(f)M*)>o - p\1 

(b 1 (t)b\(t)b 2 (t)b\(t)} = l- P \l-2p 1 i 1 
PIT = (b\(t)bt(t)) , 

pll = P 1 £ = -(b 1 {t)b\{t)) . 



pll 
pll 



(6) 



At the thermodynamic limit N = oo, we get the following analytical formulations, 
(&I(t)&i(t)>o = ^- f (l-4,alP 2 k s\n 2 A k t)dk, 

i r n 

(b\{t)b 2 {t)) = — / cosfc(l - 4a 2 k (3ls\n 2 A k t)dk, 



(bi{t)b 2 {t)) 



1 

2^ 



sin£;a fc /3 fc (2sin 2 A fc i(l - 2f3\ 



is\n2A k t)dk. 



(7) 



We use concurrence [19|] to measure the bipartite entanglement. When the density matrix 



is the "X" type, it is given by 



obtain t 
discord 



20| 2max[0, \p%\ - y^pE, \p\l\ - ^fp%pg}. In order to 
re dynamics of quantum correlation and classical correlation, we choose quantum 



11] to measure the quantum correlation (it equals to the quantum correlation 21] 



in the case of two qubits [14J, |2l|, |22j). It is given by the discrepancy of two kinds of quantum 
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FIG. 1: Evolution of (a) entanglement, (b) quantum correlation and (c) classical correlation for 
the nearest sites, with anisotropy parameter 7 = 1.0. 

mutual information derived from their classical counterparts, D{pvz) = I{pvi)— C(pi2), where 
C(Pi2) = max {n 2 }(>S'(Pi) ~~ Ylii Q-j^iPx))- {n?} are the complete projective measurement on 
partition 2, and qj = tri^O^jpn)- p{ = (H^p^H^/qj is the state of partition 1 after partition 
2 recording the outcome j. 

The initial state of the spin 1 and spin 2 reads 

(\ o\ 






P12 



(8) 



\0 0/ 

so there are neither classical correlation nor quantum correlation between 1 and 2. As the 
magnetic field quenched (A changes suddenly from zero to a finite value), the initial state will 
evolve under the XY exchange interaction. In Figure 1(a), (b) and (C), evolution of entan- 
glement, classical correlation and quantum correlation have been shown with the anisotropy 
parameter 7 = 1.0. The behavior of entanglement is complex. There is a definite boundary 




FIG. 2: (Color online) The first local maximum of classical correlation for different values of 
anisotropy parameter. The solid black, dash red, dot blue, short dash dark cyan line denote 
7 = 1.0,0.9,0.8,0.7 respectively. 
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FIG. 3: The evolution of (a) entanglement, (b) quantum correlation and (c) classical correlation 
for the next-nearest neighbor sites with anisotropy parameter 7 = 1.0. 
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around A& ~ 0.86, when A < A& entanglement is present in a long time interval, but when 
A > Xb entanglement disappears fast and absents for a longtime. However, the quantum 



3, 



correlation is nonzero almost all the time. So it is clear that it is due to the dissonance 
another kind of quantum correlation, the quantum correlation could be nonzero without en- 
tanglement. From the evolution of both entanglement and quantum correlation we could see 
that when the magnetic field is large (v), there will be entanglement all the time. While for a 
little magnetic field, the dissonance, will be the dominant component of quantum correlation 
and evolves with time. The classical correlation behaves much the same like the quantum 
correlation, as the magnetic field decreases classical correlation oscillates much faster, it can 
be understood as the magnetic field decreases the flipping interaction dominates over the 
on site potential energy. It can be seen that for all kinds of correlation, they disappear 
gradually as the magnetic field dying out. This is due to the fact the exchange interaction 
tend to dispose the nearest sites in the sparable state (| f)i + | 4)i) <8> (| f) 2 — | 4)2) or 

(lt>i-ll>i)®(lt> 2 + U> 2 ). 

For different values of the magnetic field, entanglement, quantum correlation and classical 
correlation evolve with time, and they will obtain their first local maximum with different 
values. Particularly, we find that for the classical correlation, its first local maximum C max 
has a maximum around the critical point A c = 1. This may be seen as an indicator of the 
quantum phase transition at the critical point. In fact we find that when the anisotropy 
parameter is big enough, 7 > 0.7, C max always has a maximum around the critical point, as 
shown in Fig. 2. But as the anisotropy parameter decreases, the maximum gets away from 
the critical point. 

For the next-nearest neighbor sites, we only need to focus on sites 1 and 3. The reduced 
density matrix takes on the same shape as that of sites 1 and 2, with 3 instead of 2 for the 
diagonal elements, and the non-diagonal elements are given by 

Pll = P^ = -(&iW^W)o + 2((6i(t)^(t)) (& t 1 (t)&iW)o 

+ (6i(t)4W)o + |(&iW& 2 (t))o| 2 ), 
p\l = Pli* = -(fei(t)& 3 (t))o + 2((6 1 (t)6 3 (t))o(& t i(t)&i(t))o 

+ 2Re((b 1 (t)bl(t)) )(b 1 (t)b 2 (t)) , (9) 

where the transverse symmetry has been used, and Re means the real part of the number. 
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The correspondent correlation functions are 

(b{(t)b 3 (t)) = ^ J cos2k(l - 4a 2 k p 2 k s\n 2 A k t)dk, 

(bi(t)b 3 (t)) = ^ J sm2ka k f3 k {2sm 2 k k t{l-2(5l) 

-ism2K k t)dk. (10) 

Here for sites 1 and 3 the initial state is also represented by Eq. (7), the dynamics of 
entanglement is obviously different from that of nearest neighbor sites, for large anisotropy 
parameter the entanglement does not appear until after a finite time interval, as shown 
in Fig. 3(a). This is a little counterintuitive, because entanglement between sites 1 (2) 
and 2 (3) comes into being almost instantly due to their direct interaction. Therefore, it is 
reasonable for the entanglement between sites 1 and 3 increases from zero once sites 1 (2) and 
2(3) are entangled. At the same time, there is a small band of A depending on the anisotropy 
parameter, in which no entanglement will be generated. What is more, for both classical 
correlation and quantum correlation, far away from the critical point, they become larger 
as the magnetic field diminishes, as shown in Figures 3(b), 3(c). This is a little unbelievable 
at the first glance, and is obviously different from the case of the nearest neighbor sites. At 
this case we do not find that the dynamics of entanglement, quantum discord and classical 
correlation could be used as indicator of the quantum phase transition at the critical point 
like the nearest neighbor sites. 

In this paper, we have studied the dynamics of entanglement, quantum correlation and 
classical correlation with one dimensional transverse XY model for both the nearest and next- 
nearest neighbor sites. We found that for the nearest neighbor sites, when the magnetic field 
is strong, all three kinds of correlation are larger than those with a weak field. While for the 
next-nearest neighbor sites, the dynamics is obviously different. For the nearest neighbor 
sites, two kinds of dynamics are shown with different values of the magnetic field, while 
for the next-nearest one, there is a narrow band of A, in which there is no entanglement 
generated. What is more, we found that for the nearest neighbor sites, the first maximum 
C max of classical correlation peaks around the critical point for large anisotropy parameter, 
it may be used as an indicator of quantum phase transition. 
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